Spectral Line Width Broadening from Pair Fluctuation in a Frozen Rydberg Gas 
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Spectral line width broadening in Rydberg gases, a phenomenon previously attributed to the 
many-body effect, was observed experimentally almost a decade ago. The observed line width 
was typically 80-100 times larger than the average interaction strength predicted from a binary 
interaction. The interpretation of such a phenomenon is usually based on the so-called diffusion 
model, where the line width broadening mostly originates from the diffusion of excitations. In this 
paper, we present a model calculation to show that diffusion is not the main mechanism to the 
line width broadening. We find that the rare pair fluctuation at small separation is the dominant 
factor contributing to this broadening. Our results give a width of about 20-30 times larger than 
the average interaction strength. More importantly, by turning off the diffusion process, we do not 
observe order of magnitude change in the spectral line width. 



PACS numbers: 32.80.Ee,34.20.Cf 



Rydberg gases have attracted renewed interest in re- 
cent years due to the unprecedented advancement in laser 
cooling and trapping [J 0, i, fl H, & 0, 1, @- Rydberg 
atoms, possessing a large dipole moment and long life- 
time, can interact with each other coherently for rela- 
tively long times, which make it a potential candidate 
for quantum information processing fiol [ll[ • There have 
been many experiments exploring the quantum many- 
body effects in Rydberg gases, e.g. spectral line broad- 
ening [H, 0, [1], number correlation @, [l3| and col- 
lective excitation Q. Many experimental results can be 
understood from the well-known dipole blockade effect: 
when two Rydberg atoms are close enough, the dipo- 
lar interaction will shift them out of resonance with the 
external driving laser, thus double excitation is greatly 
suppressed. 

In this paper, we are interested in the unusual line 
width broadening which was observed in experiments 
[H, 0. To be specific, we will consider the following 
two cases. (I), the main process is \np) + \np) «-> 
\ns) + \ {n + l)s) for experiment where the principal 
quantum number n is 23 and the maximal gas density 
is around 10 10 cm~ 3 . np, ns, and (n + l)s are abbrevi- 
ated as p, s, and s' , respectively. (II), the main process 
is |(n + l)s) + \n's) «-> \np) + \ (n' + l)p) for experiment 
0, where the principal quantum numbers n and n' are 
24 and 33, respectively, and the maximal gas density is 
around 10 9 cm -3 for each of the s state, (n + l)s, n's, 
np, and (n' + l)p are abbreviated as s, s' , p, and p' , 
respectively. For both cases, they express the creation 
process, e.g., in case (I), one atom makes a downward 
transition from the Rydberg state \p) to \s), and the 
other atom makes an upward transition from \s') to \p), 
that is to say, creating ss' from a pp pair. The detuning 
between \pp) and \ss') is controlled by a static electric 
field and the transition is allowed with dominant dipole 
moments /i ps and /i ps ' . Here u a p denotes the transition 
dipole moment between states \a) and |/3). Similar no- 
tations will apply to case (II). In addition to the above 
creation processes, there also exist the exchange process, 
e.g. \p) + \s) «-» |s) + \p). Different from creation process, 



the exchange process is always resonant and it describes 
the hopping of excitation in the whole gas. For this rea- 
son, we will also call it a diffusion process. 

A rough estimate from binary interaction will give a 
line width of the order of the average interaction strength 
V\ = u ps u ps rn or V2 — H ps u p i s rn for the two cases respec- 
tively, with n the average density of the gas. However, 
the experimentally observed line widths are typically 
~ lOOVi for the first case, and are ~ 8OV2 for the sec- 
ond case. Previous explanations are based on a diffusion 
model, where resonant processes like \p) + |s) <-> \s) + \p) 
form a diffusion band and play a dominant role in the 
broadening. In this model, even at large detuning, there 
are still some pairs of atoms close enough to perform 
the creation process. Because the hopping of excitation 
can happen in the whole gas, the diffusion evacuates the 
excitations so that each pair of close atoms can react sev- 
eral times, analogous to autocatalytic processes in chem- 
istry. As pointed out in Ref. the band formed by 
the diffusion of \ss') is coupled to the state \ss') and the 
existence of this band broadens this population transfer 
0, , which shows that the broadening should be a re- 
sult of the many-body effect. However, we believe that 
it is mostly a two-body effect arising from the density 
fluctuation, as we will show below by simple theoretical 
reasoning and numerical simulations. 

To better understand the diffusion process, we consider 
a toy model for the process \s) + \p) \p) + \s). The 
Hamiltonian is 



H 



Tk) 



(1) 



which is purely a diffusion process under the dipolar in- 
teraction 



Vdip(^) = Cd 



1-3 cos 2 



(2) 



r = ri — r 2 , cos 9 = z ■ f, and Cd = uj p . The energy scale 
is chosen to be V = nl p n, the inverse of which sets up the 
time scale. To ease our discussion, we consider only one s 
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atom and a bunch of p atoms with zero magnetic moment, 
i.e., magnetic quantum number m = 0. In Fig. [TJ we 
show the probability of finding an s state on the initial 
s atom as function of time, i.e., P — \(ip(t)\ijj(t = 0))| 2 
where \ip(t)) is the many-body wave function at time t. 
We can see that the probability decays smoothly and 
saturates to a finite value. At t ~ 0.2, there is about 
50% chance that the s state has drifted away. So the 
characteristic time of diffusion is on the order of 1/10. In 
addition, we show the histogram of eigen-energy in the 
inset of Fig. [T] by directly diagonalizing the Hamiltonian 
(H|). We find that the width of this diffusion band is 
roughly 5, corresponding to 1/t. This poses questions 
on the original explanation of the band diffusion model. 
Since for large detuning, e.g. A = 40 Vi, in order to make 
a non-negligible transfer from pp to ss' , the interaction 
strength between them should be of the same order as 
A. In this case, pp and ss' should be split by an amount 
of ~ A. That is to say, the manifold of pp and ss' will 
be in a large detuning to the band. Thus, the state ss' 
with a large detuning from the band will not decay into 
it. Therefore, the explanation of broadening from this 
diffusion band model is questionable. 




t(1/V) 

FIG. 1: The probability of finding s state on the initial s 
atom as function of time. The inset shows the histogram 
of the eigen-energy. The simulation is carried out with 256 
atoms and over 1000 spatial configurations. Time is in unit 
of 1/V (see text). 

We investigate the line width problem of case (I) and 
(II) by direct numerical simulations. To do this, we ran- 
domly put several atoms in a cubic box and assume each 
atomic state \p) or \p') has no magnetic moment. The 
numerical schemes are as follows. We first calculate the 
full dynamics with a given detuning for each spatial con- 
figuration of atoms with a fixed evolution time. We then 
average over spatial configurations to obtain the excita- 
tion probability as function of detuning from which we 
can extract the line width. In order to minimize the finite 
size effect, we adopt wrap boundary condition. 

We first focus on case (I) and discuss case (II) subse- 



quently. The energy and time scale are the same as those 
in the toy model. In case (I), each atom can be in the 
states p, s, or s' . The Hamiltonian is found to be 

H = Y}y^e'^ t \p ] p k )(s ] s' k \ + V; k \p ] s k ){s. ] p k \ 
jk 

+ V^\ Pj s' k )(s' jPk \} + h.c. 
describing the following processes 



p + p 
p + s 



P 



s + s , 
s+P, 
s'+P, 



(3) 
(4) 
(5) 



where process ((3|) is not always resonant and its de- 
tuning A is controlled by an electric field, while pro- 
cesses dU and ([5]) are always resonant. Vj k , VL, and 
Vj' k all take the form of V<a p (fj — f k ) with corresponding 
Cd = ^s P ^sp',P-%, and iij. p , for processes ([3]), Q and (JSJ), 
respectively. Initially all of the atoms are in the state p, 
and they evolve under the dipolar interaction for a fixed 
time T. We are interested in the yield of s atom (/ s ) as 
function of detuning, from which we can extract the line 
width. 

We perform calculation with up to N = 10 atoms and 
average over the initially random atom positions many 
times. Under wrap boundary condition, our results al- 
ready show convergent behavior for N = 8 atoms. Ex- 
trapolation to N = oo will give about 15% difference to 
that of N = 10. As we will see later, this difference is 
not crucial to our conclusion as we are interested in the 
order of magnitude difference. Our results are shown as 
solid lines in Fig. O For the marked solid line, we in- 
clude all three processes, while for the unmarked solid 
line, we only consider the process ([3]), i.e., we effectively 
turn off the diffusion process. The extracted line widths 
are about 35 and 25 for the two cases, respectively. This 
shows that the process © already gives a width of order 
of magnitude larger than the average interaction strength 
and the diffusion process further broadens the line width 
by roughly 50%. So the unusual line width broadening 
mostly comes from the process ©• 

In real experiments, the density is not uniform. There- 
fore, we need to take into account the density profile of 
the atomic cloud. We assume it as a Gaussian form with a 
width a, so the density is written as n(r) = 2\/2ne~ r ^ a 
where n is the average density. Accordingly, the excita- 
tion fraction with a Gaussian convolution (/ s q) is given 
by 



/* 



s,G 



J f s (A/V x V/V{r))n(r)r 2 dr 
J n(r)r 2 dr 

Jf s {A/V x e" 2 /" 2 /2V2)e- r2 ^ 2 r 2 dr 



(6) 



Here we assume the density varies slowly on the length 
scale we considered. The results for N = 10 atoms are 
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shown as the dash-dotted lines in Fig [2l The marked 
dash-dotted line is the result including all three processes. 
While the unmarked dash-dotted line is the result only 
including process ([3]). The extracted line width with dif- 
fusion is about 30, which is a few times smaller than the 
results in Ref. [lj. The one extracted from the calcula- 
tion without diffusion is about 20, again demonstrating 
that diffusion is not the main mechanism in the broad- 
ening. We also note that the curves are not perfectly 
symmetric around A = due to anisotropy of the dipo- 
lar interaction. However, this difference is too small to 
be detected under current experimental conditions. 
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FIG. 2: (Color online) f s and f s ,a as function of A for N = 10 
atoms averaged over 1000 spatial configurations. Results 
with marked (unmarked) lines include (exclude) diffusion. 
Solid and dash dotted lines are for the case of homogeneous 
and Gaussian convolution, respectively. Parameters used: 
T = 3.4, a = 500, ju SJ) = 1.02, (i^ = 0.98. The results 
are insensitive to a. A is in unit of Vi (see text). 

We further consider the motional effect on line width 
broadening. To do this, we give a constant speed v s for 
each atom but with random direction. For v s = 0.05, the 
atom moves 0.1 (average distance) at the end of simu- 
lation,which is still in the so-called "frozen gas" regime. 
Our results give an additional broadening of about 20%, 
so the motional effect is not important as expected. 

Following similar procedures, we investigate case (II). 
Now each atom can be in the states s, s', p, or p' . The 
Hamiltonian is found to be 

H = Yl{ V jke- iAt \PjPk)(sA\+VjkM(sm\ 

jk 

+ V-' k \p' j s' k ){sjp' k \] + h.c. 
describing the following processes 



s + s 
p + s 
p' + s' 



p + p', 
s+p, 
s' + p' 



(7) 
(8) 
(9) 



with dipolar interaction. Our numerical results are 
shown in Fig. [3] In this case, the width with (with- 
out) diffusion is found to be about 20 (10). Therefore, 
the diffusion does not cause order of magnitude change 
in the line width broadening. 
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FIG. 3: (Color online)/ p and f Pt c as function of A for a to- 
tal N = 20 atoms averaged over 1000 spatial configurations. 
Results with marked (unmarked) lines include (exclude) dif- 
fusion. Solid and dash dotted lines are for the case of homo- 
geneous and Gaussian convolution, respectively. Parameters 
used: T = 0.36, a — 500, /J, sp = 2,/Vj>' — 0.5. The results 
are insensitive to a. The initial population of s and s' are the 
same. A is in unit of V2 (see text). 

So what causes the broadening? It is nothing unusual 
but the rare pair fluctuation at small distances. To see 
this point, we calculate P(|A|), which is the probability 
distribution of nearest neighboring atoms with the ab- 
solute interaction strength not larger than the absolute 
detuning. P(|A|) can be found with the help of 



P(|A|) = 



|A| dP(\V\ < |A'|) 
d|A'| 



d|A'| 



(10) 



In a homogeneous system, according to the Erlang dis- 
tribution [l5j], the nearest pair distribution is cx e - 4vr3 / 3 
with unit density. Therefore, for an isotropic interaction 
Viso = 1/r 3 with Cd = 1 (e.g. V = 1), we have 



dP(|Vi so | < |A|) _ 4tt 



d|A| 



3IAI 



(11) 



While for dipolar interaction Vdi p with a = 1 (e.g. V 
1), we have 



dP(\V dip \ < |A|) _ Air 



d\A\ 



3IAI 



efo|l - 3a; 2 |e 3 i A 



1 1— 3a; 2 



(12) 

As |A| — > +00, the asymptotic behaviors are 
(4tt/3)|A|- 2 and (16tt/9\/3)| A|" 2 for isotropic and dipo- 
lar interaction, respectively. As |A| — ► + , they approach 



4 



and \/3/(47r) for isotropic and dipolar interaction, re- 
spectively. The remarkable difference at |A| — > + is a 
signature of the dipolar interaction. 
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FIG. 4: (Color online) P(| A |) as function of |A|. Blue dashed 
lines are for isotropic interaction and red solid lines are for 
anisotropic interaction. A is in unit of V (see text). 

For |A| = 40, the probability of nearest atom pairs 
that has an interaction energy larger than |A| is about 
10%, i.e. those atom pairs will have non-negligible contri- 
bution to the dynamics. The calculated f s of case (I) at 
this detuning is about 8%, close to the estimated value. 
Therefore, the rare pair fluctuation is the main cause of 
the line width broadening and the diffusion of excitation 
further increases this broadening by roughly 50%. 

For case (II), we can also compute the line width (w) 
for different ratio of s and s', i.e. io as function of v = 
(m — + n^) where 711(712) is the density of s(s') 

atoms. Our numerical results are shown in Fig. [SJ The 
largest error in our simulation still comes from the finite 
atom effect which has been discussed for Figs. 2 and 3. 
This error does not change much as we vary the ratio of 
s and s' atoms. Other errors are negligible. We find that 
w increases as v increases and saturates at v = ±1. The 
increasing behavior of w with v is due to the imbalanced 
hopping of s and s' atoms (n sp = 4/Vp/). The increased 



ratio in s atoms will thus show a stronger diffusion effect. 
However, this increase in the line width is again not an 
order of magnitude change. 

To conclude, we have reexamined the important role 
that pair fluctuations play in the spectral line width 
broadening of a frozen Rydberg gas. From direct numer- 
ical simulations, we find that density fluctuations con- 
tribute to a width of roughly 20 — 30 times of average 
interaction strength. In addition, by turning off the diffu- 
sion process, we did not find order of magnitude change in 
the line width. Therefore, the large line width is primar- 
ily due to density fluctuations, and the diffusion process 
is not overwhelmingly dominant as previously suggested. 
However, our numerical results are only in qualitatively 
agreement with experimental results. The even larger 




FIG. 5: Line width w as function of population ratio v of s 
and s' atoms. The same parameters are used as in Fig. [3] 
except for different ratio of s and s' atoms. 



line width observed in both cases and the double peak 
structure as observed in (I) cannot be explained from 
current calculations which encourages further investiga- 
tions. 
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